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The intrinsic magnetic moment (IMM) and intrinsic angular momentum (lAM) of a chiral 
superconductor with p-wave symmetry on a two-dimensional square lattice are discussed on 
the basis of the Bogoliubov-de Gennes equation. The the IMM and lAM are shown to be 
on the order of /isY and HN, respectively, N being the total number of particles, without 
an extra factor (Tc/Tp)'^ (7 = 1, 2), and parallel to the pair angular momentum. They arise 
from the current in the surface layer with a width on the order of the coherence length ^0i the 
size of Cooper pairs. However, in a single-band model, they are considerably canceled by the 
contribution from the Meissner surface current in a layer with the width of the penetration 
depth A, making it difficult to observe them experimentally. In the case of multi-band metals 
with both electron-like and hole-like bands, however, considerable cancellation still occurs for 
the IMM but not for the lAM, making it possible to observe the lAM selectively because the 
effect of the Meissner current becomes less important. As an example of a multi-band metal, 
the case of the spin-triplet chiral superconductor Sr2Ru04 is discussed and experiments for 
observing the lAM are proposed. 


1. Introduction 

1.1 Issues concerning intrinsic angular momentum in ^He A-phase 

The size of the intrinsic angular momentum (lAM) Lin has been a major issue since the 
middle of the ’70s, when the size of the lAM in the A-phase of superfluid ^He was keenly 
discussed. Namely, the issue is the value of the exponent 7 when the lAM in the ground state 
is expressed as Lin = {Alh/2){Tc/Tp)'^ in the case where the Lvector is uniform in space. 
This problem was first addressed by Anderson and Morel in their seminal paper^) discussing 
a possible state of superfluidity of liquid ^He, who stated that 7 = 1. This estimation is based 
on the idea that the lAM is sustained by a pair condensate with angular momentum h so that 
the lAM is proportional to the pair amplitude (oc Tc/Tp), as is clearly elucidated by Leggett in 
a review article on superfluid ^He.^^ After the discovery of the superfluidity of ^He, a bunch of 
theoretical works were performed on the size of the lAM by microscopic calculations based on 
the k-space approach, which predicted 7 = This result seems to have been interpreted 
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such that the lAM is formed by two particles inside the thin surface layer in a fe-space near 
the Fermi surface whose width is on the order of the superfluid gap A, implying that the 
number of particles forming the lAM is of the order of A x (A/Flp)^ or N x (Tc/Tp)^. 

On the other hand, in 1976, Ishikawa proposed an alternative idea®) that the lAM in the 
A-phase is Ljn = {Nh/2) in the ground state if we adopt the wave function proposed by BCS, 
whose property was discussed by Ambegaokar in the review article appearing in the textbook 
on superconductivity edited by Parks.®) The point is that the orbital-part of the many particle 
ground state ^(ri, r 2 , rs, r 4 , • • • , ryT-i, ttv) [N being an even number) is given by the product 
of the wave function of the Cooper pair (p{r, r') which is the eigenstate of the relative angular 
momentum = h concerning (r — r'). Namely, 

T(ri,r2,r3,r4, • • • ,rN-i,rN) = A [(/^(ri, r2)¥5(r3, r4) • • • (/?(rAr_i, fat)] , (1 ) 


where A indicates the anti-symmetrization with respect to all the spatial coordinates, and the 
spin coordinates are discarded as irrelevant. Indeed, Ambegaokar showed explicitly that this 
form of the wave function can be transformed into the wave function proposed by BCS if the 
Fourier component of the wave function of the Cooper pair (/?(r, r') is well defined with respect 
to the relative coordinate (r — r'). The issue at that stage was whether the last condition is 
satisfied or not. Indeed, the Fourier component of — r'), is given by 


— ) 
Uk 


( 2 ) 


where the coefficient Uk expressing the coherence of the BCS state vanishes outside the thin 
region around the Fermi surface of the width of cut-off Ec so long as the weak-coupling theory 
is applied. However, the weak-coupling theory might discard a crucial effect on the coherence 
between the region near the Fermi surface and the core region of the Fermi sphere. Namely, 
the issue is to what extent the coherence of the Cooper pair condensate is maintained into 
the core of the Fermi sphere in a real situation. This seems to be a fundamental question 
regarding how to understand the Fermi superfluid state, including superconductivity. 

It is instructive to remember the discussions of Bogoliubov et a/.^®) and Anderson and 
Morel,^) who estimated the effect of the Coulomb repulsion on the Cooper pair formation. 
They solved the gap equation with the phonon-mediated attractive interaction near the Fermi 
surface with a width on the order of the Debye energy wd and the Coulomb repulsive inter¬ 
action acting in the whole fe-space, and they obtained a superconducting gap that is finite 
not only in the thin surface around the Fermi surface but also in the whole fc-space including 
the core of the Fermi sphere. This was crucial to understanding the systematic deviation of 
the index a for the isotope effect from the canonical value a = 1/2 Therefore, it is a real 
effect that the coherence of the Cooper pair extends to the core of the Fermi sphere. In this 
sense, a difficulty posed for the idea of Ishikawa was safely avoided.Adopting the wave 
function in Eq. (1), McClure and Takagi showed^^) that the result Lin = Nh/2 holds more 
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generally within the manifold of the ground state with an axially symmetric configuration 
of the ^-vector, such as the Mermin-Ho^^^ and Anderson-Toulouse^®^ textures. Namely, the 
state with a uniform configuration of the ^-vector is adiabatically continued from such axially 
symmetric states. 

A related issue concerned the structure of the supercurrent js in the A-phase of superfluid 
^He. By the symmetry argument, js is expressed as 

js = Ps'^s + (5(V X 1) + —V X , (3) 

where ps is the superfluid density tensor and the tensor C has the following form in general: 


Cij = C5ij-CoUr (4) 

The issue concerned the form of C together with the size of Ljn. According to microscopic 
calculations based on the k-space representation, Cij is given as Cij = C{T){6ij — 2£iij) and 
C{T = 0) = hn/4, where n is the number density, and Ljn ~ husiTc/Tp)'^^’^'^') or Lm = 0,^i 
where Ug is the superfluid number density. On the other hand, the supercurrent density at 
T = 0 is given by 


1 

Js = PsVs + 2^ ^ 



( 5 ) 


according to a modified wave function based on Eq. (1), which is generalized so as to take into 
account a gradual variation of the center of mass coordinate in This supercurrent 

is in a form parallel to that appearing in the electromagnetism in materials,in which the 
electric current jM induced by a variation of the magnetization density M is given by 

jM = — V X M. (6) 

ho 

Indeed, if one remembers the relation M = (eh/2m)L, L being the orbital angular momentum, 
and that js corresponds to (m/e)jM, the second term in Eq. (5) is precisely Eq. (6). On the 
other hand, Mermin and Muzikar^®) claimed that the calculation in Ref. 12 misses a subtle 
singularity of the wave function (p(r, r') when it is calculated in the k-space representation, 
leading to the expression for the supercurrent 

js = pvs -h ^/i(Vns) X e. (7) 

This form was also derived by Nagai on the basis of kinetic theory.The physical meaning 
of the difference between Eqs. (5) and (7) has been discussed by Volovik from a general 
viewpoint.In any case, it is crucial that both expressions, Eqs. (5) and (7), include the 
term proportional to (Vug) x i that makes an essential contribution to the lAM through the 
surface current when i is parallel to the surface where the number density n vanishes abruptly. 

This fact leads to the following physical picture of the lAM in a uniform configuration of 
the A vector. In the region at a distance from the system boundary of more than the coherence 
length ,^0 of the Cooper pairs, the relative motions of the Cooper pairs around the Avector 
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cancel with each other, resulting in no contribution to the lAM. On the other hand, in a thin 
region with a width on the order of near the system boundary, its cancellation is incomplete, 
giving rise to the surface current that is the main origin of the I AM. This physical picture 
is exactly the same as the picture for an electric current induced by the (classical) orbital 
magnetism in materials. Precisely speaking, it is not self-evident that Eq. (5) or (7) cannot 
be applied near a system boundary where the number density changes abruptly, because they 
were derived on the assumption that the spatial variations of physical quantities are gradual 
compared with the length scale characterizing the physics, i.e., the coherence length ^o- 

Recently, this problem has been revived in the context of the topological effect associated 
with the surface state of a chiral superfluid or superconductivity. According to detailed cal¬ 
culations that take into account the microscopic structure of the surface state with spatial 
variation on the order of the result derived from the supercurrent, given by Eq. (5) or (7), 
is essentially correct.On the other hand, these calculations show that the lAM at finite 
temperatures is given by 


Ain(T) 


Ni\T)h 


( 8 ) 


where Ns (T) is the superfluid number parallel to the £-vector.This result cannot be easily 
understood intuitively, but seems to be much more involved than a naive physical picture. 
In particular, it is difficult to find the physical reason why appears in a two-dimensional 
system when the ^-vector is perpendicular to the two-dimensional plane. Indeed, it should 
be compared with the result obtained by a calculation using a cylindrical representation for 
one-particle states, in which L\^iT) is given by 


Ain(T) 


Ns^{T)n 

2 


(9) 


where N^{T) is the superfluid number perpendicular to the £-vector. 

Regarding an experiment for detecting the lAM in the ^He-A phase, NMR measurement 
of the structural change in the Mermin-Ho texture in a rotating cryostat was proposed by 
Takagi.^^^ This experiment has been performed at the Institute for Solid State Physics (ISSP) 
of the University of Tokyo and suggested that the size of the lAM is on the order of Nsh/2?^'^ 


1.2 Issues concerning intrinsic magnetic moment in Sr 2 Ru 04 

In the past decade, the problem concerning the lAM has been revived as that of the intrin¬ 
sic magnetic moment (IMM) in the spin-triplet chiral superconductor Sr 2 Ru 04 ,^®^ in which 
the orbital part of the superconducting gap has been identified as = A(sin, /i: 3 ; 0 -|-isin kyo), 
a being the lattice constant of the two-dimensional lattice. This is also supported by a muon- 
spin-resonance (/iSR) measurement showing the breaking of time-reversal symmetry.^*^) Infor¬ 
mation on the fc-representation of A^ is obtained from the temperature dependence of the 
specific heat (under a magnetic field) and theoretical investigations that suggest the impor- 
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tance of short-range ferromagnetic correlations among quasiparticles.If the lAM Ljn is 
on the order of N^h/2 and the gyromagnetic ratio is given by {—el2m), with e(> 0 ) being 
the elementary charge, as in the classical case, the IMM density is estimated as 

( 10 ) 


ns m 

^ '''■band 


where Us = Ns/V, = 4 ttx 10 “^ H-m“^ is the magnetic permeability of vacuum, fiB = eli/ 2 m 
is the Bohr magneton, and ^ is the harmonic average of the band mass of electrons over 
the occupied state in the Brillouin zone discussed in Appendix A and should be distinguished 
from the effective mass averaged over the Fermi surface discussed in Appendix B. . Then, the 
magnetic flux density without the external magnetic held H is given by Min because the 
relation B = M -\- holds by dehnition.^®^ The electron number density n of the 7 -band, 
which is the electron-like band, in Sr 2 Ru 04 , is estimated as 

1 


n = 


abc 


( 11 ) 


where a = b = 3.9 x 10“^® m and c = (12.7/2) x 10“^® m are the lengths of an edge of the 
primitive cell of Sr 2 Ru 04 along the a, b and c directions, respectively.The magnetization 
density Mjn is given by the relation 

e 


Min— occ — ^nfiQ- ^ , 


( 12 ) 

^band 

where ^ ~ 2.9 m is the effective mass of the 7 -band of Sr 2 Ru 04 .^^^ Therefore, the intrinsic 
magnetic hux density Bin is estimated as 

10-30 


Bin = — - 


m 


nhc oyiOcc 

"^band 


X 5.8T 


-2.1 X 10”^ T = -2.1 X 10^ G. 


(13) 


This value is larger than the “observed” lower critical held Rex’® = 5.0 x 10“^ T for Sr 2 Ru 04 .^^) 
Therefore, at hrst sight, this intrinsic magnetic hux density Bin will not be completely screened 
by the Meissner current. However, since Sr 2 Ru 04 has two other bands, a hole-like a-band and 
an electron-like /3-band, considerable cancellation in the IMM is expected among the electron¬ 
like /3- and the 7 -bands and the hole-like a-band, as discussed in Sect. 5. 

Another issue is whether the lAM is also cancelled by the orbital angular momentum aris¬ 
ing from the Meissner current moving through a thin surface with the width of the penetration 
depth A near the boundary of the system. If this cancellation is incomplete, the remnant lAM 
may be detected by the Richardson-Einstein-de Haas effect.^®’ 


1.3 Purpose and organization of the present paper 

In the “rotationally symmetric” system, in which the angular momentum is a conserved 
quantity, the result for the lAM, Lin = {Nh/2)i, is almost self-evident from the viewpoint 
of the BEC-BCS crossover or of the adiabatic continuation. In the BEG limit, each diatomic 
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molecule has angular momentum M as shown in Fig. 1(a). Therefore, the lAM is given by 
Lin = where the number of diatomic molecules is N/2. So long as the pairing 

interaction has rotationally symmetry, the value of Lin, which is a conserved quantity, should 
not change even if the pairing interaction is weakened to approach the BCS limit in which 
“molecules” overlap each other as shown in Fig. 1(b). A subtlety is that the gap amplitude 
A of a Cooper pair has a very weak singularity. Namely, cPA/d^'^ has a discontinuity as the 
chemical potential /i passes through the bottom of the quasiparticle band.^®^ 

On the other hand, in a lattice system in which the rotation symmetry is broken, the 
problem is not so trivial. Therefore, it is necessary to explicitly investigate the problem for 
a specific lattice model. One of the purposes of this paper is to investigate how the results 
for the lAM of the ®He-A phase in three-dimensional free space are modified in the case 
of a chiral superconductor on a two-dimensional square lattice, which simulates Sr 2 Ru 04 . 
Another purpose is to investigate, using this two-dimensional model, to what extent the IMM 
is screened by the Meissner effect and the I AM is lessened by the Meissner current. On the 
basis of our results, it is discussed how the IMM and lAM are observed in Sr 2 Ru 04 . 



Fig. 1. (a) Picture of strong-coupling (BEC) limit, where diatomic molecules with angular momentum 

M are in the state of Bose-Einstein condensation, (b) Picture of weak coupling (BCS) limit, where 
the angular momentum is distributed on the Cooper pair condensation. 
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The organization of the present paper is as follows. In Sect. 2, we introduce the model on 
a square lattice with an attractive interaction between nearest-neighbor sites and discuss a 
formulation for explicit calculations. In Sect. 3, the results for the lAM and IMM are shown. 
In Sect. 4, the effect of the Meissner current on the size of the IMM and lAM is discussed. 
Finally, in Sect. 5, we propose how to observe the IMM and lAM of the two dimensional chiral 
superconductor Sr 2 Ru 04 . 


2. Chiral Superconductor on Square Lattice 

2.1 Model Hamiltonian 

In order to study the problem of the lAM and IMM of a chiral superconductor on a 
two-dimensional lattice, a model of Sr 2 Ru 04 , we start with the following Hamiltonian: 

'^ = -l^Yl - t (14) 

i(7 {i,j)a (i,j)a 

where /r, t, and V are the chemical potential, the transfer integral between nearest-neighbor 
(n.n.) sites of the square lattice, and the attractive interaction between electrons at n.n. sites, 
respectively, and c|^ (cia) is the creation (annihilation) operator of an electron at the i-th site 
with spin component a (=t or ),). The symbol {i,j) indicates that the summation is taken 
over the n.n. sites. Here, we consider the spin-triplet pairing with 5^ = 0 and introduce a 
superconducting gap in the spin-triplet manifold as 

V 

+ cqCj-i-), (15) 


where (• • •) means the average by the mean-field Hamiltonian given as 


Tdmi 


ia 


t Y 4aCja + I + 4A4 + 

(*d> ^ 



(16) 


Here the gap A^- depends on lattice sites i and j in general, and its dependence is determined 
self-consistently by solving the (lattice version of the) Bogoliubov-de Gennes equation together 
with Eq. (15).^®^ The gap Aij is odd with respect to the interchange of i j: 




(17) 


which results in the odd-parity pairing. Note that, in the case of a uniform system without a 
boundary, the most stable gap among those given by Eq. (15) is expressed in a wave-vector 
representation as^^^ 


Ak = A (sin feajO-|-i sin kyo), 

where a is the lattice constant. 


(18) 


2.2 Orbital magnetization and angular momentum in band picture 

In order to take into account the effect of the magnetic held B, we adopt the following 
way of giving the Peierls phase the transfer integral tij between electrons at the i-th and j-th 
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sites 


iijiB) = tij exp 


h 


r 

Jn 


A(r) • dv 


(19) 


where A is the vector potential giving the magnetic field as B = V x A and the contour 
integral is performed along the line connecting the two sites. Then, the band-energy part 
^band of the Hamiltonian is expressed as 

"Hband = - ^ iijiB) c\^Cja. (20) 

(*uV 

Adopting the gauge of the vector potential such that A = (B x r)/2, "Hband is reduced to the 
following form: 

^band = -t Ei [exp Olfyia) + exp (-i|f 2 /*o) 

E* [exp (i||yio) + exp {-^^Via) 

Ei [exp (i||xia) c{.^ + exp (-i|fxia) 

+t E* [exp (i|f x^a) + exp (-i|f 4.,’ (^1) 

where the spin coordinate is abbreviated for concise presentation and the i-th position of the 
lattice is designated by the two-dimensional representation {ix-,iy)- In deriving Eq. (21), the 
contour integral Eq. (19) along the x-direction has been approximated by the trapezoidal rule 
as 


rxi 

J Xi—a 


Axix,yi)dx 


and 


l‘Xi-\-a , Q, \ 

J Ax{x,yi)dx Ax\^Xi +-,yij a. 


Xi- -,yi}a 


( 22 ) 


(23) 


A similar approximation is adopted for the integral along the y-direction. 

Then, the operator of the magnetization M of the system (parallel to B, z-component) is 
given, in the limit H —)• 0, as 


Ito 


dn 
dB ) 


/ <9^band \ 

V dB J 


ta 


B=o 


i-e)'^iri X pi)^, 


(24) 


where the “momentum” operator pj at the i-th site is defined by 

Pxi 2 a i^(ix+l,iy)a — ^(ij;,iy)ai^iix+^jy)<^ ^{ix — i,iy)<x) 

Pyi 2 a i^{ix ,iy + l)o' — ^(ix,iy)cri^iix,iy + ^)o' ^{ix,iy — i)cr') 

If we introduce the band mass rub at the zone boundary, say at the T-point, as 




1 


(25) 


(26) 


2 mb ’ 

the magnetization operator Mz [Eq. (24)] is reduced to a band version of the conventional 
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i 

This implies that the definition of the orbital angular momentum L, 

L = ^(rj X Pi), (28) 

i 

is a valid and natural one. The definition of mb [Eq. (26)] corresponds to the free-electron-like 
dispersion of tight-binding dispersion around the T-point, {kx,ky) = (0,0). Namely, 

Cfc = —2t(cos k^a + cos kya) ~ —At -|- ta^{k^ + ky) -!-•••. (29) 

Using this dispersion, m(]“^ in Eqs. (10), (12), and (13) is estimated as m(]“^ = (7r^/4)mb in 
the half-filled case, as shown in Appendix A. 

3. Results for I AM and IMM 


An explicit form of the Bogoliubov-de Gennes equation for the mean-field Hamiltonian 
[Eq. (15)1 with the superconducting gap of jSz = 0 [Eq. (18)] is given by^°^ 

eui = —/i Ui — t Uj + ^ipji (30) 

UA) 

evi = fjLVi + tvj + '^ ^ijUj , (31) 

(hi) 

where {j, i) means that the summation is taken over the nearest-neighbor sites. An actual 
calculation is performed as follows. Hereafter, we mainly focus our discussion on the half- 
filled case, unless otherwise stated. Equations (30) and (31) are diagonalized by means of a 
unitary transformation 14 to give the mean-field Hamiltonian 

Hrai = ^ ^ (32) 

m=l m=l 

where Al is the number of lattice sites, 0 < ei < ei • • • < e^v^, and the fermion operators 7 
describing the quasiparticles are related to the electron operators c by 

[^If ’ ’ ’ ’ *"AfLt’ ’ ’ ’ ’ ''AfLtl 

= [7iV> • • • > 7]vLt’ 7U> • • • , (33) 


Substituting Eq. (33) into Eq. (15), we obtain the self-consistent equation for the gap Ajj as 
Y 


V 


m=l 

JVl 


E l(^)j+Ni^,m+NL(^)i,m+NL ~ (^)i+Ni^,m+NL (^)j,m+NL] /(^ m)j (34) 


m=l 
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where U depends on Ajj and Sm {'m = 1, - ■ ■ , and f{x) is the Fermi distribution function 
f{x) = {e^ + 1 ). 

We have solved Eqs. (15), (16), and (32)-(34) self-consistently using the numerical diag- 
onalization method and obtained the gap and the energy levels Sm {'m = I,-- - ,A^l)- 
Numerical calculations have been performed for lattice sizes of up to Ap = 30 x 30 squares 
with both open and periodic boundary conditions. Throughout the present paper, the phase 
of the superconducting gap A^ is chosen as shown in Fig. 2, while Aj [i = 1-4) are determined 
self-consistently. 



Fig. 2. Phase of superconducting gap A^- with i-th site chosen as a center. 


First of all, the IMM [Eq. (27)] and lAM [Eq. (28)] are shown to be zero for the periodic 
boundary condition for which Ai = A 2 = A 3 = A 4 , resulting in the recovery of the chiral 
gap given by Eq. (18) in the k-space representation. This is because the orbital currents due 
to the relative motion of the Cooper pairs cancel with each other. On the other hand, if we 
use the open boundary condition, we obtain finite values of the IMM and lAM because the 
cancellation of the relative motion of the Cooper pairs is incomplete near the boundary of the 
system within the coherence length ^0 of the Cooper pairs, as in the case of the classical theory 
of magnetic materials where the current jm given by Eq. (6) exists only in a thin surface layer 
of the system if the magnetization M is uniform in the bulk of the system. 

Figure 3 shows the phase diagram in the T/t-V/t plane for the system size Ap = 9x9 
in the half-filled case. Three different superconducting states exist, which we call SI, SII, and 
Sill, as shown in the figure. The SI phase appears around the phase boundary between the 
SII phase and the normal phase in the intermediate-coupling region Vjt < 4.5. On the other 
hand, the Sill phase appears not only around the phase boundary between the SII phase 
and the normal phase in the strong-coupling region Vjt > 4.5, but also as a low-temperature 
phase in the strong-coupling region Vjt > 8.5. 

Figure 4 shows the pattern of current in the SI phase for the parameter set, Vjt = 3.25 
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Fig. 3. Phase diagram in T/t-V/t plane. The system size is taken as A^l = 9x 9. Three supercon¬ 
ducting phases, SI, SII, and Sill, exist together with the normal phase N. 


and T/t = 0.564, shown by a closed sqnare in Fig. 3. The system size is taken as Alp = 9x9. 
The phases of the snpercondncting gap Ajj shown in Fig. 2 are given by Ai = Ai, A 2 = —iA 2 , 
A 3 = —A 3 , and A 4 = iA 4 , where Ai ~ A 4 are real and have the same sign. Note that the 
sign of Ai ~ A4 is determined so as to satisfy Eq. (17). In this phase, the current is induced 
around the center of the bulk and its pattern can be seen as a deformed vortex pair lattice 
as shown by symbols ® and 0 . The physical reason why this pattern is realized is unclear for 
the moment. 

Figure 5 shows the pattern of current in the SII phase for the parameter set, V/t = 8.00 
and T/t = 1.482, shown by a closed circle in Fig. 3. The system size is taken as Al = 9x 9. This 
phase is the bulk phase in the intermediate-coupling region as shown in Fig. 3. The phases of 
the superconducting gap Ajj shown in Fig. 2 are given by Ai ~ A 4 , which are real and have 
the same sign. The pattern of the current distribution is that expected physically. Namely, 
the current is induced near the boundary of the system owing to incomplete cancellation of 
the relative angular momentum of the Cooper pairs. The reason why the current exists at the 
center of the system is that the system size Ap = 9 x 9 is comparable to the extent of the 
Cooper pair ^*(r) = 7 r^(r), ^(T) being the coherence length at finite T. 

If we follow the BCS theory for an s-wave weak-coupling superconductor, the coherence 
length in the ground state (at T = 0) is given by .^0 = hup/vrA. Then, the ratio of .^0 to the 
lattice constant a is estimated as 


^0 27 Ap r. 1 1 

— — - ~-~ () 1 1 - 

a air'^Tr vr^ Tc Tc 


(35) 


where we have used the BCS relation A = nTc/j, where log 7 is the Euler constant 
C = 0.57721 • • • , and assumed the free dispersion for the quasiparticles and that the Fermi 
momentum is given by pp ~ hir/a. As shown later in Fig. 8 , Tc/t ~ 1.8 for the attractive 
interaction V/t = 8.00. The Fermi energy Ep is estimated as Ep ~ 4t in the case of half-filling. 
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Fig. 4. Pattern of current in the SI phase for V/t = 3.25 and T/t = 0.564 as indicated by a closed 
square in Fig. 3. The system size is taken as =9 x 9. The current at each site is shown by 
an arrow whose length represents the relative size of the lattice momentum dehned by Eq. (25). 
Symbols 0 and © indicate localized vortices with opposite circulations. 


Therefore, with the use of Eq. (35), ^o/o is estimated as io/a ~ 0.24, giving the estimated 
extent of the Cooper pair = tt^q as ^*/a ~ 0.75. On the other hand, the temperature 
T/t = 1.482 for Fig. 5 is about 90% of the transition temperature Tc, i.e., — T)/Tc ~ 0.1, 

as seen in Fig. 3. Then, with the use of the correlation length ^(T) ~ 0.74^o\/7c/(Tc —T) 
at a finite temperature the extent of the Cooper pair at T/t = 1.482 is estimated as 

^*{T)/a ~ 0.75 X 0.74 x VTo ~ 1.76, which is not negligible compared with the system size 

= 9 X 9. 

Figure 6 shows the pattern of current in the SII phase for the same parameter set as 
above, V/t = 8.00 and T/t = 1.482, but for a much larger system size, = 30 x 30. We can 
see that the current is essentially confined near the system boundary with width ~ 1.76a. 
The current along and near the boundary of the system can be regarded as a lattice version 
of the surface current of the A-phase of ^He, which is a topological superfluid characterized 
by the topological class D defined in Ref. 43. 

Figure 7 shows the pattern of current in the Sill phase for the parameter set, V/t = 9.00 
and T/t = 0.644, shown by a closed triangle in Fig. 3. The system size is taken as Ap = 9x9. 
The phases of the superconducting gap A.jj shown in Fig. 2 are given by Ai ~ A 4 which 
are real and have the same sign. In this phase, the induced current forms concentric layers 
of flows, with adjacent layers having opposite signs. The physical reason why this pattern is 
realized is unclear for the moment. 
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Fig. 5. Pattern of current in the phase SII for Vjt = 8.00 and T/t = 1.482 as indicated by a closed 
circle in Fig. 3. The system size is taken as A^l = 9x 9. Current at each site is shown by an arrow 
whose length represents the relative size of lattice momentum defined by Eq. (25). 


In Fig. 8 , mz = which is the IMM per site divided by fiQ, and Tc/t, the 

superconducting transition temperature, are shown as functions ofV/t, which is the strength 
of the attractive interaction, at T/t = 0.01 for the system size Afp = 13 x 13 at half-filling 
{fi/t = 0). Note that, in order to avoid the effect of the boundary, Tc/t is calculated with 
the periodic boundary condition. It is noteworthy that does not decrease even though Tc 
decreases and that a negative correlation exists between niz and Tc/t except in the SI phase 
where the coherence length becomes comparable to the system size and the superconducting 
state is greatly suppressed over the whole system. This suggests that the lAM, connected to 
the IMM by Eqs. (27) and (28), is on the order of hN/2 without the extra factor (Tc/Tf)'’' 
(7 = 1 or 2 ). 

Indeed, according to Eqs. (27) and (28), the lAM, Lin, is expressed in terms of HN and 
mz/auB as follows: 


Lin — 


2mb Mz 


mz 

-hNi^- 


, . , (36) 

(-e) Ho m anB 

where a = tma?/f^ and we have used the definition of the Bohr magneton hb = eh/2m. Since 

it is easily derived that am]c/m = 1/2 if we use the definition of a given above and Eq. (26), 

the lAM at half-filling (i.e., Ap = A) is given by 

m . 

(37) 


2 an-B 


Considering the size of mz/otHB in the SII phase shown in Eig. 8 , we can see that the lAM in 
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Fig. 6. Pattern of current in the SII phase for V/t = 8.00 and T/t = 1.482. The system size is taken 
as iVn = 30 X 30. The current at each site is shown by an arrow whose length represents the 
relative size of the lattice momentum defined by Eq. (25). 

the SII phase is given by Ljn ~ hN/2. This implies that 7 = 0 for the exponent of the extra 
factor (Tc/Tp)'^, verifying the validity of the original theory by Ishikawa for the lAM in the 
^He-A phase.®) 

Figure 9 shows the system-size dependence of for the case of V/t = 3.0 and 

T/t = 0.01 at half-filling, i.e., fi/t = 0, up to A^p = 19 x 19. We can see that the system-size 
scaling works reasonably well, giving m^/a/xe — 0.302 in the limit A^p oo. 

Figure 10 shows how IMM/^q per site, niz/afiB, changes depending on the strength of the 
attractive interaction V/t, together with the behavior of A/t, with A = Ai = A 2 = A 3 = A 4 
dehned in Fig. 2. The system size is taken as Ap = 30 x 30, the maximum size adopted in 
the present paper. The reason why approaches zero at A/f ~ 1 can be understood as 

follows. At V/t ~ 1, the transition temperature is Tc/t ~ 0.08 as shown in Fig. 8 . Therefore, 
according to ^o/« given by Eq. (35) and the extent of the Cooper pair ^*/a for 

Tc/t ~ 0.08 is estimated as 

r 4t , , 

— = — ~ vr X 0.11 X -~ 17. (38 

a a 0.08f 

This implies that the superconducting order for the system size Ap = 30 x 30 is almost 
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Fig. 7. Pattern of current in the Sill phase for Vjt = 9.00 and T/t = 0.644 as indicated by a closed 
triangle in Fig. 3. The system size is taken as IVl =9x 9. The current at each site is shown by 
an arrow whose length represents the relative size of the lattice momentum defined by Eq. (25). 


m^lafiB Tcjt 



Fig. 8. IMM ///0 per site niz as a function of the attractive interaction V/t for the system size TVl = 

13 X 13 at T/t = 0.01 in the case of half-filling, i.e., /i/t = 0. SI, Sll, and Slll indicate the phases 
shown in Fig. 3. The dimensionless parameter a is defined by a = tma^/I t?. 


destroyed by the effect of the boundary of the system. 

We can see that mx/a/XB decreases as V/t increases. The V/t dependence of in 

the region V/t > 10 can be fitted by const./(Id/t), although we do not show this explicitly. In 
this region, the pattern of the current corresponds to the lattice of vortices and antivortices 
with domain walls as shown in Fig. 11 for the system size A^l = 15 x 15 with V/t = 50 and 
T/t = 0. Namely, the tendency that the magnetic moments of the vortex and antivortex cancel 
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mzlajjL-Q 



Fig. 9. System-size scaling for IMM//ro per site ruz for V/t = 3.0 and T/t = 0.01 at half-filling, i.e., 
= 0. The dimensionless parameter a is defined by a = tma? /h^. 



Fig. 10. IMM//io per site rUz and magnitude of the superconducting gap A as functions of the 
attractive interaction V/t for the system size Al = 30 x 30 at T/t = 0 at half-hlling, i.e., ^/t = 0. 
The dimensionless parameter a is defined by a = tma^/h^. The dashed curve is that expected in 
the limit Al —?► oo. 


with each other becomes prominent in this region. However, this extremely strong coupling 
region does not seem to be realized in actual systems. 

The dashed curve in Fig. 10 is a smooth extrapolation to V/t —^ 0 where mz/a^B is 
expected to be niz/a^B = d/vr^, as argued below. Indeed, according to Eq. (10) and Eq. (A-4), 
which is valid at half-filling, and the definition of ce, a = tma?/h^, mz/a^,B is transformed as 
follows: 

ruz 1 m /TB 14m 
a/iB 2 m)(“ ^ a/XB 2a vr^ mb 
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Fig. 11. Pattern of current in the Sill phase for V/t = 50 and T/t = Q. The system size is taken as 
IVl = 15 X 15. The current at each site is shown by an arrow whose length represents the relative 
size of the lattice momentum defined by Eq. (25). Symbols (8) and 0 indicate localized vortices 
with opposite circulations. 


4 

= — = 0.40528---~ 0.41, (39) 

where, in deriving the last equality, we have used the relation amb/m = 1/2, which is derived 
from Eq. (26). 

Figure 12 shows the filling (uc) dependence of IMM//ro, for the attractive 

interactions V/t = 2.0 and 3.0 at T/t = 0. The reason why rriz /approaches zero at 
ric — 0.15 in the case oiVjt = 3.0 can be understood by the effect of competition between 
the system size and the extent of the Cooper pairs as in the case in Fig. 10. The extent of the 
Cooper pairs is estimated by Eq. (35) but with a replaced by a/^/n/ because the average 
distance between electrons increases in inverse proportion to the square root of the filling ric- 
Then, the following relation holds: 


— ~7rx0.11 


T ^ 

The transition temperature is calculated as — 0.071t. Therefore, for V/t = 3.0 and 
the filling ric = 0.15, which gives Ep ~ 0.887t, is estimated as 

e 


~ 11 , 


(41) 


which is approximately half of the size of the system of 30a. This explains why rriz becomes 
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zero at approximately Uc = 0.15. 

The dashed curve in Fig. 12 is a smooth extrapolation to ric 0, where mz/ajiBUc is 
expected to be equal to 1. Indeed, in the dilute limit [ric 0), the effect of the lattice fades 
away so that the result in free space should be recovered. The lAM Lin in tbe free space is 
expected to be given by Lin = fiN/2. On the other hand, extending Eq. (37), Lin is given by 


Ln = -5iV 


m. 


where = N/N^. Therefore, in the limit ric 


(y.fi'Qfic 

0 , rriz/auBnc is expected to approach 1. 


(42) 


mz/a/^Bric 



Fig. 12. (Color online) IMM/^xq per particle rriz/auBnc as a function of the filling ric for the attractive 
interactions, V/t = 2.0 and 3.0. The system size is taken as = 30 x 30. The dimensionless 
parameter a is defined by a = The dashed curve is that expected in the limit ric ^ 0. 
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4. Effect of Meissner Current on lAM 


In this section, we discuss the effect of the Meissner current on the IMM and lAM in¬ 
duced in a chiral superconductor as discussed in the previous section. Here, we estimate the 
orbital angular momentum due to the Meissner current flowing in a thin surface layer within 
penetration depth A. First of all, let us consider the situation shown in Fig. 13(a), where the 
flat boundary surface of the superconductor is the xy-plane and the magnetic field Bx{z) is 
parallel to the x-direction and decreases in the superconductor (z > 0). Then, by the Ampere 
law, the current density jy{z) is given by 

1 dBx{z) 


jyiz) = 


(43) 


fiQ dz 

Therefore, the Meissner current Jm (per unit length), flowing in the ^-direction in a surface 
layer in the x-direction, is given by 

1 dBx{z) _ Bq 
/^o ’ 


poo 

Jm= dz- 

Jo f^o dz 


(44) 


where Bq = Bx{0). The corresponding mass current Pm is estimated as 


tpc 

Pm —- Jm, 


(45) 


where is the band mass averaged near the Fermi level. The reasons why the mass 

appears are that the actual current is caused by a deformation of the Fermi surface, and that 
the Fermi liquid correction by “P®” due to the back flow effect does not cancel the dynamical 
mass enhancement in the case where the Galilean invariance is broken in lattice systems such 
as Sr 2 Ru 04 .‘^'^^‘^®i Then, in the situation shown in Fig. 13(b), the orbital angular momentum 
Lm due to this Pm is given by 

Pm = 271R X R X Pm- 


With the use of Eqs. (44) and (45), Pm is reduced to 

Pm = 2'kRJ X 


,2 .. "ibid 


P, 


X - 


e V Po. 

If we equate Pin = in Eq. (12) to Pq at the surface. Pm is finally given by 

_ n(7rp2)n _ Nh 

— - ^ 


m 


band 


9 

^ "^band 


(47) 


(48) 


The hrst equality of Eq. (48) is also valid for columnar systems with a general cross-section 
shape if the factor vrP^ is replaced by the cross-section area S because the angular momentum 
is related to the areal velocity. In this sense, the expression for Pm, Eq. (48), is valid for 
columnar systems with an arbitrary shape of the cross section. This result implies that Pm 
cannot cancel the I AM P^ = Nh/2 in general as far as / "^band’ which is usually the 

case including the case of free space, as for liquid ^He, if the many-body effect near the Fermi 
surface is taken into account. Namely, the total angular momentum P in the ground state. 
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given by 


L — + Lm 


Nh 

~Y 



^hLd \ 

^occ / 


(49) 


is not vanishing but is on the order of 0{Nh) and has the opposite sign to the i vector 
because ^ in correlated systems in general. We note, however, that the total 

angular momentum L, given by Eq. (49), vanishes in the hypothetical system with a free 
electron dispersion because "iband/”^band = 1 in such a case. 


(a) 




Fig. 13. (Color online) (a) Magnetic field near the boundary between a superconductor and a vacuum, 
(b) Distribution of the magnetic field (indicated by arrows) in a cylindrical superconductor sample. 


5. Multi-band Effect 

The discussions so far has been for a single-band model. However, multiple bands exist 
in general. In particular, Sr 2 Ru 04 , which is a promising candidate for exhibiting the lAM 
and IMM, has three bands, one hole-like band (a) and two electron-like bands (/3 and 7 ) (see 
Fig. 14). Therefore, the IMMs of the electron-like and hole-like bands partially cancel with 
each other because they have opposite signs, reflecting the opposite signs of the gyromagnetic 
ratio. The characters of the bands are summarized in Table According to Eq. (10), the 
IMM is inversely proportional to the band mass. Therefore, the relative size of the IMM is 
also inversely proportional to the band mass. The relative value of the IMM is given in the 
last row in Table I for the case that the lAMs of each band have the same magnitude. Then, 
total IMM density, M^°^, is approximately given by 


A-C ^ ( "n + ^ ^ 


n. 


1 


1 


I 


(50) 


if Eq. (10) is valid for the three bands shown in Table I. The factor in parenthesis in Eq. (50) 
gives a small value of — 0 . 110 ; thus, the magnitude of is one order smaller than that 
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Table I. Band-dependent mband/™ according to Ref. 47 and the lAM and IMM ratios. 


Band 

Character 

a 

Hole 

/5 

Electron 

7 

Electron 

"Iband/m 

1.1 

2.2 

2.9 

lAM ratio 

1 

1 

1 

IMM ratio 

1 

-0.5 

-0.38 


expected in the single-band model: 

Ti 

~ X 0.110. (51) 

Therefore, MT* induces the magnetic field B\^ ~ 2.3 x 10“^ T = 23 G, which is smaller than 
the lower critical field ~ 5.0 x 10“^ T and is easily screened out by the Meissner effect, 
in contrast to the case of the single-band model as discussed in Sect 1.2. 



Fig. 14. (Color online) Fermi surfaces of Sr 2 Ru 04 . The a-band is hole-like, and the /3- and y-bands 
are electron-like. 


On the other hand, this cancellation does not occur in the total lAM, is given 

by the sum of the lAM of three bands, which are all expected to be on the order of 0{Nsh/2). 
Therefore, the effect of the Meissner current on the total lAM is much less important than 
in the case of the single-band lAM given by Eq. (49). Thus, the total lAM remains to be 
technically unscreened by the Meissner current. This multiband effect is crucial to distinguish 
the lAM from the IMM. 
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6. How to Observe IMM and lAM in Sr 2 Ru 04 

In this section we discuss how to detect the IMM and lAM by carrying out experiments 
that are possible to perform on Sr 2 Ru 04 . 

As discussed in the previous sections, it will be difficult to detect the magnetic field induced 
by the IMM owing to the Meissner screening in the bulk sample. However, it may be possible 
to detect it for a small sample with size on the order of the penetration depth A r\j 1300 A. 
On the other hand, the current pattern near the boundary of a small sample will be rather 
complicated in the sense that the directions of the Meissner current and the current inducing 
the IMM are opposite and the ranges of both currents are different, leading to considerable 
cancellation between the two currents up to the penetration depth A. Therefore, in order to 
detect the IMM, the experiments should be carefully designed. 

Another possibility exists of observing the IMM by performing /rSR experiment, which can 
detect the IMM induced around the stopping site of in principle. This is because attracts 
electrons from adjacent Ru sites and acts as a nonmagnetic impurity that destroys the chiral 
superconducting order at sites surrounding /r'*', resulting in a local circulating current around 
the site because the cancellation of the chiral current of Cooper pairs becomes incomplete 
there. If the effective impurity potential is sufficiently strong to completely suppress the chiral 
order at surrounding sites,the magnetic field induced by this circulating current will be on 
the order of Rin [Eq. (13)], which can be easily shown by standard calculations in classical 
electrodynamics based on the expressions for the surface current given by Eqs. (5) and (7). 
A crucial point here is that at the /i'*' site is free from the Meissner screening effect. 
On the other hand, if the effective impurity potential is not so strong, the induced magnetic 
field at the site is decreased considerably depending on the strength of the potential. 
Indeed, the magnetic field observed by /rSR, — 0.5 G, is rather small compared with 

Hin — 2.3 X 10“^ T = 23G given by Eq. (51). This fact can be understood by assuming that 
the effective impurity potential from is only moderate. However, this should be verified by 
an explicit model calculation, which is left for a future study. 

On the other hand, the I AM of each band gives additive contributions to the total I AM on 
the order of 0{Nsh) which can be probed by the so-called Richardson-Einstein-de Haas effect, 
which has been used to detect the macroscopic spin angular momentum in the ferromagnetic 
state.Since the size of the lAM is on the same order as the spin angular momentum 
of ferromagnetic compounds, it is expected that the lAM can be observed by this effect in 
practice, although a sufficiently low temperature will be required. 

7. Conclusion 

On the basis of the tight-binding model with the nearest-neighbor attraction on a square 
lattice, the IMM and lAM have been calculated by solving the Bogoliubov-de Gennes equation. 
It turned out that, in the ground state, the IMM per site divided by fiQ is on the order 
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of //B and the lAM is on the order of ~ HN. In particular, in the dilute limit, the result 
strongly indicates that the I AM approaches Lin = hN/2, which is the value first predicted by 
Ishikawa.®) The lAM and IMM are induced by the surface current flowing in a thin layer with 
a width on the order of the coherence length ^o- 

It has been shown that thus created lAM is partially screened by the Meissner current if 
the particles have an electric charge. The extent of the cancellation depends on the ratio of 
the effective mass of quasiparticles near the Fermi level to that of the harmonic average 

"^band occupied State: if they were equal, the cancellation would become perfect. This 

is in marked contrast to the IMM, which is almost completely screened by the Meissner effect 
if the spontaneous magnetic field created by the IMM is smaller than the lower critical field 
Hci, as in the case of Sr 2 Ru 04 . 

On the other hand, it turned out that the multiband effect is important, as in the case 
of Sr 2 Ru 04 . Namely, a certain amount of cancellation of the IMM occurs among particle and 
hole bands because they have charges with different signs, while such a cancellation does not 
occur for the I AM. 

An interpretation of the spontaneous magnetic held observed in Sr 2 Ru 04 by /iSR, and a 
possible means of probing the lAM have also been proposed. 
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Appendix A: Harmonic Average of Band Mass of Electrons on Square Lattice 

In this appendix, we derive the harmonic average of the band mass of electrons with 
dispersion, Eq. (29). The inverse mass tensor m~^ of a band electron is given by 

( d^ek \ 

dkl dk^dk,, 


-1 


m = 


1 

¥ 


d^ek 


dkxdky 




dky 




\ dkydkx 

Substituting the dispersion of an electron, Eq. (29), into this expression, it is easily seen that 
the matrix in Eq. (AT) is diagonal. Indeed, its explicit form is 



cos kytt 


The averaging of ihxx over the occupied states at half-filling is performed as follows: 


{^XX ) 


2t¥ 


1 


rnia 


i(7r/a)2 Jo 


dk.. 


rw/a ky 2taJ 4 

dkx cos kxtt = —TTw X —rr 


'0 


h? 


TT^ 


(A- 2 ) 


(A-3) 


where the area of integration with respect to kx and ky is restricted to the part surrounded 
by the dashed line (the Eermi surface) of the first quadrant in Eig. AT. The averaging of khyy 
is performed in a similar way, giving the same value. Then, using the definition of the band 
mass (near the T-point) mo [Eq. (26)], the harmonic average of the band mass over the 

occupied states is given by 


moJnd = 


1 

— = —mo- 


{jdxx ) 


(AT) 
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Fig. AT. Brillouin zone of electrons on the square lattice. Dashed lines indicate the Fermi surface at 
half-filling and dotted lines indicate the order of two-dimensional integration. 


24/27 



















J. Phys. Soc. Jpn. 


Full Paper 


Appendix B: Effective Mass Averaged over the Fermi Surface 

In this appendix, we derive the effective mass averaged over the Fermi surface of electrons 
with dispersion, Eq. (29). The inverse mass tensor of this band electron is given by Eq. 
(A-2). Therefore, the average of m~]. over the Fermi surface is calculated as follows: 


i^xx ) 


2ta 


2 J dkcos kxa6{cos kxa + cos kya + ///2t) 


dk(5(cos kxtt + cos kya + /i/2t) 


1 2ta2 

2 ~W 


dk (cos kxa + cos kya)6{cos kxa + cos kya + fJ-/2t) 
J (ik(5(cos kxa + cos kya + fi/2t) 


1 fia^ 

'2'W' 


(B-1) 


The expression for {iTiyy) is also given by Eq. (B-1). Therefore, the effective mass averaged 


over the Fermi surface is given by 


"^band — (B-2) 

It is remarked that diverges toward half-filling, i.e., /r —)• 0, which is consistent with the 

existence of the van Hove singularity in the density of states at // = 0. On the other hand, 
in the dilute limit, i.e., fi —>■ —4t, approaches k?/2ta?, which is the same as the band 

mass mb near the F-point defined by Eq. (26). This guarantees the validity of the definition 
in Eq. (B-2). In the almost-filled case, ^ ~ 4t, approaches —h^/2ta? = —mb, the hole 

band mass near the top of the band. 
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